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Controlling the time evolution of the population of two states in cavity quantum electrodynamics (QED) 
is necessary by tunning the modified Rabi frequency in which the extra classical effect of electromagnetic 
field is taken into account. The theoretical explanation underlying the perturbation of potential on spatial 
regime of bloch sphere is by the use of (Bagrov, Baldiotti, Gitman, and Shamshutdinova) BBGS-Darboux 
transformations^ on the electromagnetic field potential in one-dimensional stationary Dirac model in which 
the Pauli matrices are the central parameters for controlling the collapse and revival of the Rabi oscillations. It 
is shown that by choosing o\ in the transformation generates the parabolic potential causing the total collapse 
of oscillations; while {02,(73} yield the harmonic oscillator potentials ensuring the coherence of qubits. 
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I. INTRODUCTION 

Recent research on quantum theory scrutinizes the rec- 
onciliation between quantum and classical physics, and 
it attempts to take the advantages from the classical 
for constructing the new theories^. For instance, the 
quantum information between sender and receiver can be 
teleported by involving classical information transmission 
generated by the joint measurement of EPR pair 6 . In the 
modern research of experimental quantum information, 
the notions of Coulomb and photon blockades represent 
a recent resurgence of interest in the influence of classical 
in quantum physics^^. Here, we show that concurrent 
condition of the electromagnetic field behaving classically 
and quantumly in cavity QED provides the conditions for 
steering the dynamics of collaps and revival sequences of 
Rabi oscillations. 

Rabi oscillations are the important features in the field 
of experimental quantum information, since the oscil- 
lations represent the coherence and decoherence of the 
qubits in various type of physical schemes of quantum 
computer 9 i 10 i 26 i 38 . 

The famous Rabi oscillations are derived from the 
quantity associated as the atomic inversion, W(t), de- 
fined as the variation in the excited- and ground-state 
populations 2 ^ 

W(t) = P e (t) - P g (t). (1) 

When the electromagnetic field is classical, the oscilla- 
tions in atomic inversion have constant amplitude. If the 
field is quantum, the oscillations consist of a sequence of 
collapses and revivals. The object of the control of the 
time evolution of the population of two states in cavity 
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is to formulate the Hamiltonian H, which realizes the 
population quantity of excited- and ground-state 2 ^. 

In the recent years, there have been several efforts to 
control the collapses and revival o 12 ! 30 ' 32 ! 37 , in order to 
damp ou t 13 ' 22 ' 33 or to ensure the stability 3 ^ of the oscil- 
lations by considering the trade-off of physical constraints 
in the related quantum computer scheme. Nevertheless, 
the general theory of how to control the atomic inversion 
is poorly comprehended. 

On the theoretical side, Darboux transformations con- 
tributions into the wide range of physics fields have been 
widely known over centuries 2 ^. Different with the other 
transforms such as the Fourier and Laplace transforms, 
the action of Darboux transform does not change the 
domain of the potential and state of the system, but it 
creates the new potential and state of the system. The 
output of Darboux transform is a set of new potential 
and state of the physical system. For instance, in nonlin- 
ear Schrddinger equation, Darboux transform generates 
the new solutions of the equation 2 ^. The Darboux trans- 
form is very useful to describe the evolution of a physical 
system due to the change of its Hamiltonian potential as 
found in wide range of nonlinear physics. However, its 
contribution to the field of quantum computation is still 
poorly understood. 

Recently, there are several efforts in order to open the 
mystery of Darboux transformations into quantum com- 
putation. Samsonov et al., followed by Hussin et alJ£, 
propose a novel technique to apply Darboux transforma- 
tions into Jaynes-Cumming Hamiltonian 3 -^. Bagrov et 
al. coined a Darboux transformation which does not vi- 
olate the two-level system equation structure^. These 
processes can be done by using the method of inter- 
twining on the exactly solvable Dirac potentials which 
are equivalent to the Darboux transformations between 
Schrodinger potentials 3 -. 

The application of Darboux transformations to one- 
dimensional Dirac equation is performed by intertwining 
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the Hamiltonian in Dirac equation 



Lh = hit, (2) 



where L is intertwining operator and ho (hi) is the old 
(new) Dirac Hamiltonian. In the sense of perturbation 
theory, ho (hi) is unperturbed (perturbed) Dirac Hamil- 
tonian, where hi = h + V(a, /?), in which {a, /?} con- 
tributing into the perturbation terms. Recently, it has 
been shown that numerous problems in quantum me- 
chanics can be solved under this technique^. 

To date, the closest theoretical result to the experi- 
mental was only successfully performed by Samsonov et 
al. by showing the contribution of pseudosymmetry po- 
tential in controlling the qubits dynamics of supercon- 
ducting circuits based on Joscphson tunnel junction^. 
In this paper, we would like to show that it is possible to 
control Rabi oscillations in cavity QED by considering 
the presence of extra classical electromagnetic field. It 
emerges the perturbation effect on the system and can 
be represented by one fold Darboux transformations on 
one-dimensional Dirac equation. The results suggest that 
there is a correlation of Dirac potentials and the atomic 
inversion. 

The paper is organized as follows : in Subsection (jl Ap , 
we introduce the method proposed by Bagrov et al., then 
we generalized the methods : Pauli matrix in operator 
B is not constrained in 03, but the use of {a 1,02} is 
also considered. Moreover, the potential is not only re- 
stricted in (pseudo)-scalar form, but also we allow the 
potential in the form of vector where the Pauli matri- 
ces act as the potential basis and the function as the 
potential coefficient. Subsection pB[) contains the in- 
troduction of the mathematical model of atom-field in- 
teraction in cavity QED. Here, the Rabi frequency and 
oscillations of fully-quantum mechanical model are mod- 
ified due to the presence of the extra classical effect of 
electromagnetic field. Section (jlT| provides the elucida- 
tion of the implementation of BBGS-Darboux transfor- 
mations to semiclassical Rabi model. The outcome is 
to obtain the parameters in BBGS-Darboux transforma- 
tions : a(t), (3(t), and |V(i)|. The section covers two 
types of procedures: first, reformulation of the equation 
from Schrodingcr equation into Dirac equation. Here, the 
representation of data is changed, but we do not control 
anything. Second, the action of Darboux transformation 
on the Dirac equation. In this step, we perturb the sys- 
tem to control the potential and state of the system. In 
Section (JTTTJ) , it is shown that the choice of Pauli matri- 
ces in BBGS-Darboux transformations relates to the type 
of one-dimensional stationary Dirac potential of electro- 
magnetic field : {(Ji} generates the parabolic potential 
(so-called Yukawa-Rabi oscillations) and {(72,(73} yields 
the harmonic oscillator potential. 



A. BBGS-Darboux transformations 

In this Subsection, we introduce (Bagrov, Bal- 
diotti, Gitman, and Shamshutdinova) BBGS-Darboux 
transformations^ which has the notation T>(oi). The 
symbol oi is the Pauli matrix in element of intertwin- 
ing operator L = + B, A = er =(J ?), and B = 

oi(t) + (/(«)- A (t)W 

The action of T>((Ji) on a set of potential and 
state of a Hamiltonian is defined by T>(ai)[N] {V, \&} 
= {V[iV], ^[iV]} transforming the old Hamiltonian 
h(V)if> = £0* -> h(V[l])^[l] = ei*[l] -> ... -> 
h(V[N})^[N] = e N i$[N}. The potential is a vector in 
three dimensional Euclidean Bloch sphere using Pauli 
matrices as the basis and it can be extended into n- 
dimensions^i. The transformations affect the potential 
from V(t) = EjvUYfi® V[m = <r(j)'7j(t), 
where /((t) = f 3 (t) + Afj(t) and a(j) E {a u a 2 ,a 3 }. 
Based on this formulation, the magnitude of poten- 
tial is transformed \V(t)\ = -> \V t [N](t)\ = 
E ■ |/j(£)|. The effect of the transformations to the state 
is ^ — > = ty[N\. In order to obtain the coefficients 
cti(t) and Pi(t), intertwining operation, £ft. id = h ncw L, 
between the old and the new Hamiltonian is needed. For 
iV-fold BBGS-Darboux transformations, Lho = hiL —> 
Lhi = h^L — > ... — > L/in-1 = h^L. 

After the above definitions, we can now express the 
following theorem. 

Theorem 1. Let {V, \P} represents a physical sys- 
tem and D(<7i)[N] is a BBGS-Darboux transformation 
operator. The open-loop control mechanism can be con- 
structed, where the BBGS-Darboux transformation oper- 
ator can be assumed as a controller, {V, VP} as the initial 
system, {V[N],^[N]} as the final system, and the eigen- 
values £jv(Ci:) are the output variables. 

Proof. Let {V, is the initial condition of a physical 
system, in which the initial eigenvalues £0 are belong to 
the system. 

The action of BBGS-Darboux transformation on this 
system is V(a l )[N] {V, *} = {V[N],V[N]}, in which 
the eigenvalues £jv(cj) arc belong to the new eigenstates 
#[AT]. 

The complete proofs of the above theorem are given in 
SectionlHland Section lHTl In this paper, we only consider 
N = 1 or one fold Darboux transformations. 

B. Cavity QED model 

Following the Ref^i, the system of a two level atom 
and photons in a single cavity mode, if the electromag- 
netic field is classical, is represented by H = hwo^f + 
hn[a+b(t) + o-"&(i)t] where b(t) = b e~ lult . If the elec- 
tromagnetic field is quantum, the Hamiltonian is H = 
hujo^f + hujB^B + h^a+B + a-B^], where 6 (fit) i s anni- 
hilation (creation) operator and I7I 2 is the average pho- 
ton number. 
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It is shown that the atomic inversion of 
Jaynes-Cummings Hamiltonian is represented by 

W H'ZflihM = ^^^o^Fd - 2^(n + 

sin 2 (Q n t) 



Z^n=0 



1), k = 4, and 



where r2„ = \/il(K* + n 
A = ujq — uj is the resonance condition. In this paper, 
we use n=30, A = 2y/2, and 0=1. 

For our purpose, we modify the Rabi frequency by in- 
volving additional term representing the classical effect 
of electromagnetic held. The modified Rabi frequency is 



9.'; 



^n(n 2 +Wb + n+ 1). 



(3) 



The superscript 'm' denotes the modification. The 
equation (J3)) represents the condition that the electro- 
magnetic held concurently behaves classically and quan- 
tumly : the term VWb corresponds to classical field 
and vn+1 is the consequence of quantum effect of the 
field. This term can be easily derived from a Hamilto- 
nian driven by a classical light source as mentioned in 
Ref42. Another way, it also might be obtained by trans- 
forming Jaynes-Cummings Hamiltonian as explained in 
Refi 3 -^. Also, the modified atomic inversion is 



W(t) 

bl 



<7,tM*)l7,t) ' 
( 7 4|^(0)| 7 ,I) 



Y ^V(l - 2(n)V& t & + " + l 
^— ' III v 



n=0 



sin 2 (n™t) 



)(4) 



Darboux transformation contributes to the b'b term. 
For iV-fold Darboux transformations, it changes b'b — ► 
bi[l]^bi[\] -> ... ->• b\[N]bi[N], where subscript i corre- 
sponds to the choice of Pauli matrix in the intertwining 
operator. 



II. RABI MODEL UNDER DARBOUX 
TRANSFORMATIONS 

In this Section, we present the theoretical explanation 
underlying the perturbation effects of atom- field in cavity 
due to the presence of the classical fields. The perturba- 
tion causes the emergence of new potential in which it 
is perturbed on y-axis of Bloch sphere. To do so, it is 
accomplished through two step of procudures: first, the 
representation of the system is changed from Schrbdinger 
equation into one-dimensional Dirac equation. The mo- 
tivation of this change of representation is because one- 
dimensional Dirac equation is considered more promising 
rather than Schrodinger equation, since one-dimensional 
Dirac equation admitting vectors in their potential, in- 
stead of Schrodinger equation which its potential is in the 
scalar form only. Second, the action of Darboux trans- 
form is applied on the Dirac equation, physically meaning 
that the physical system is perturbed so that the poten- 
tial and the state of the system can be controllable. 

We rigorously evaluate BBGS-Darboux transforma- 
tions in Rabi model, the results substitute into the mod- 



ified Rabi frequency to obtain the dynamics of Rabi os- 
cillation. The evaluation of BBGS-Darboux transfor- 
mations in Rabi model involve three Pauli matrices : 
{ci f o~2, 0-3}. The outcomes of this evaluation are : first, 
the perturbation terms, i.e., the coefficients Qj(t) and 
j3i(t) which are required to obtain the transformations 
operator, B; second, the new Dirac potential of field, 
|Vi[l](t)|. Moreover, we also obtain the computational re- 
sult following the theory; to accomplish it, several meth- 
ods are needed to simplify the differential equations form 
and also to obtain the exact solutions of the equations. 
We find that Homotopy Perturbation Method (HPM)22, 
is necessary when solving the coupled-nonlinear Riccatti 
differential equations problem for T>(ai)[N = 1], while 
for D({<72, 03}) [-/V = 1], the choices produce ordinary 
differential equations only; however, one needs the The- 
orem in Ref . 27 ' 31 to simplify the equations in complex 
number form by finding their square roots, in case they 
can be easily simulated. 



A. Rabi model in one-dimensional Dirac equation 

The first procedure, i.e., the change of representation 
from Schrodinger equation into Dirac equation, is given 
in this Subsection. Rewrite Schrodinger equation into 
one-dimensional stationary Dirac equation of Rabi model 



where h = (icr z j- t +V), V 
Ve 

l , and 



= £ *- (5) 
-hn(a+b{t) - cr-tf(t)) and 



%e 



2 



£0 18 ^f 1 - We also use the following assumption b(t) 



bi[l](t) = 2Pi(t)-b(t). 



(6) 



Because bo is a constant and real, = bo, then. 

The Dirac potential in the equation §5§ can be eas- 
ily changed into V(t) = iMlbo(a x sin(u}t) — <j y cos(u)t)) 
meaning that the potential is in vector form using Pauli 
matrices as orthogonal basis. 

Although the idea of representing the physical system 
in one-dimensional Dirac equation may be better than 
Schrodinger equation, since the potential of Dirac Hamil- 
tonian is in vector form, however it has the drawback: the 
atomic excitation term in the Hamiltonian is vanished in 
the new representation. One of justification for this draw- 
back is that the new representation may be valuable to 
explain the phenomenon in which strong coupling is much 
greater than atomic excitation as observed in Reb 2 ^. 



B. The action of BBGS-Darboux transformation on Rabi 
model in one-dimensional Dirac equation 

Following the procedure of representation is the action 
of Darboux transformation on the new representation. 
Below, the procedure is elucidated deeply. 
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The one fold BBGS-Darboux transformations change 



V(t) ^^.1^1, = cj{lff 1 {t)+a{2ff 2 {t) 
i=i 

= -hn(a + b(t) -a-b\i)) 

= ihil,bo(a x sin(uit) — o~ v cos{uit)) 

-+ Vi[t] >>•.,!'/; = a(l) 1 f[(t)+a(2) 2 m 
i=i 

= -nn(a + bi[i](t) - a-b\[l](t)) 

= —iMl(a y (2Pi(t) — bocos{ujt)) + a x (bosin(u)t))) . 

It is clear that the transformations change the initial 
electromagnetic field over the y-direction on the Bloch 
sphere in the magnitude of 2/3j(t). Due to the transfor- 
mations, the potential magnitude changes 

\v(t)\ = (nm Q ) 2 -> 

\Vi[l](t)\ =4f3f(t)~4(3 l (t)b cos(ut) + bl (7) 

(One can find by the similar manner that for Jaynes- 
Cummings model \V(t)\ = (hil) 2 {n + 

From intertwining operation Lh = hiL, one can obtain 
the following equation: 

ia z B + Ml(a + bi[l](t) + a~bl[l](t)))B - hfl(Ba + b(t) + 
Ba~b ] {t))) - hil(a + b(t) + a~b\t))) = 0. (8) 

This is the master equation for the following Subsec- 
tions. 



1. 2?((7l) 



difficult to obtain exact solutions of coupled nonlinear 
equations using separately homotopy theory or pertur- 
bation theory. Most perturbation methods face problems 
with parameters: the methods admit the existence of a 
small parameter, in contrary there is no small parame- 
ter in frequent nonlinear problems. On the other side, 
most homotopy methods fail in defining the deformation 
of most coupled nonlinear differential equations into sim- 
pler equations due to their complexities. 

As pointed out by Sweilam et ai, the modified HPM is 
an effort to cover some drawbacks in conventional HPM. 
In conventional HPM, the solution is truncated in a se- 
ries in which often coincides with the Taylor series. The 
problems arise since the series have very slow convergent 
rate, therefore the solution under HPM may be inaccu- 
rate. For that reason, a new idea in this modified method 
is proposed by involving Padc approximant to enlarge the 
domain of convergence of the solution in which truncated 
in the Taylor series. 

To do so, the truncated solution in Taylor series as 
found in conventional HPM is transformed by Laplace 
transformation, followed by finding the Pade approxi- 
mant, and finished by taking the inverse Laplace trans- 
form to obtain a more accurate solution of the problem. 
We try to apply these similar procedures in order to ob- 
tain the exact solutions of the equation (9). The complete 
algorithms are given in Appendix (A). 

Starting from solving the system in the equation (9) 
by using conventional HPM, one can find that the linear 
equations in terms of p, which is an embedding parame- 
ter, of these two equations are 



P 



(«i)oM 



0, 



2fih 

-0i)o[t]=O, 



(10a) 
(10b) 



First, we consider T>{o\) to determine cxi(t), /?i(i), and 
|Vl[l](t)|. By substituting B = a x (t) + i(b(t) ~ (3i{t))ai 
into the equation ([5]). it produces the coupled nonlinear 
Riccati differential equations as following 



ai(t) 
2MI 



P 2 i(t) + Pi(t){e iut + b e~ wt ) + b Q = Q, (9a) 



x(e^ + e^6 ) + ((/3i)oW) 2 , 

(&)i[t] = ie- iwt uj(Sm-l)b 
+2(ai) [t]Hn((Pi) [t] - cos[tot}), 



(11a) 



(lib) 



-Pi(t) + iub e- iut (hn - 1) + 2 ai (t)Ml 
x(Pi(t)-coa(ut))=0, (9b) 

which are subject to the following initial conditions 

ft(0) = Wi)o = b and oi(0) = ( ai ) = ^fcffl - 
In order to obtain the exact solution of these equations, 
we follow the suggestion for using the Homotopy Pertur- 
bation Method (HPM)22. 

HPM is considered as a promising tool to solve exactly 
coupled nonlinear equations since the method success- 
fully reconciles the homotopy theory with perturbation 
theory^. Prior to their amalgamation, it was generally 



p2 . (^M = (e 



2m 

+209 1 )o[t](/3i)i[i], 



iujt i g — iu)t 



bo)Wi)i[t] 



03i) 2 [t] = 2n/*(ai)i[t](08i) o [t] - cos[ojt}) 
+2Qh(a 1 ) [t]{^) l [t}. 



(12a) 



(12b) 



Furthermore, the series solution obtained by using the 
iV-th order perturbation of HPM can be truncated as 
following 
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«l(*) S 



iujb (l-Ml) ofc0 ^ 2 u \n , A 1 ~ eiuJt )\ 



i=0 



2hn(b - 1) 



2Ml(bf ) - b ){t + i- 



(13a) 



TV 



/3 x (t) « ^(/?i)i(i) = 6o + (Ml - 1)((1 - e~^) 



(6o - 1) v w 



(13b) 



Then, the procedures involving Laplace transforms s is replaced by j. The definition of Laplace transform 
change jS?[ai(i)] = ai(s) and «5f [/?i(i)] = /3i(s) in which JSf[ai(i)] = c?i(s) can be found in the Appendix O 

These procedures yield 



/ \ \-^, \ , \ iuibn(\ — Ml) „^„„9 , n 9/ 



=o 



2Ml{b - 1) 



(1-ut)' 



(14a) 



N 



h(t) ~ £>i)<(t) = M + (tti - l)< 2 ( r J ^7 + 7T^( 

i=0 



b t 2 )) + O(t 3 ). 



(14b) 



+ wt (6 -l) l + w 2 t 2 

I 

(t) with {M N}> ^ G I s ' ^ as ^ S ^ G P * s obtaining the inverse Laplace 

{«(*), /3(t)} transforms, «5f[{ai(s), /?i(s)}] = {ai(i), to 

and M+N< 6, is applied on the equations (fT4|) to enlarge obtain the true solutions as following 
the convergence of these solutions followed by replacing 



Pade approximant, 



ai(t) « — (6 - bl)Ml(3ioj(4itu> + t 2 uj 2 - 2) (Ml - 1) + AMl 2 (b - l) 2 t(9itcj + t 2 uj 2 - 12)) + 0(t 3 ), (15) 



2i(26o-D(Bn-l) _ i(36o-2)e-'-*(Bn-l) _ i^CMW) _ _ j) + 2(& _ 1)t{b + m _ j) 

/Mi) « H a ^ 7T 2 — + 0(t 3 \l&) 

I 

The Dirac potential of this transformation can be ob- 2Ml(a. 3 (t)(j3 3 (t) — bocos(ujt)) — ibo(3 3 (t)(cos(ujt) + e~ %Mt 



taincd by substituting the equation (|16[) into the equation 
([7]). The plots of Rabi oscillations and potential versus 
time of this transformation are shown in the figure (p} . 

For the following Subsections, we find that the differ- 
ential equations are easier and similar. 



+i(P 3 (t)) 2 + -y (6o(l + e- l2u ") + b u(sin(cjt) + cos(U))) 



(17b) 



2. V(a 3 ) 



Second, we evaluate T>(a 3 ) to obtain a 3 (t), P 3 (t), and 
|V 3 [l](t)|. By replacing B = a 3 (t) + i(b(t) - /3 3 (t))a 3 into 
the equation (j8]), it yields 

I3 3 (t) + ia 3 {t) + iLub e~' luJt = 0, (17a) 



3. T>(a 2 ) 

Third, we examine T>(a 2 ) determine a 2 (t), fi 2 (t), and 
|V 2 [1](*)|. By changing B = a 2 (t) + i(b(t) - (3 2 (t))a 2 into 
the equation ([5]), it results 
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■ W(t) I 



(b) 



FIG. 1: The time evolution of the population of two 
states ( 1 (a) ) and the time evolution of Dirac potential 
( 1(b) ) in cavity quantum electrodynamics of T>(p\) for 
t = 100. 



2hn((/3 2 (t)) 2 - /? 2 (t)6 e- Mt ) + id 2 (t) 
-(iujb e- iut + $ 2 (t)) =0, 

—hQbo(iu))(sin(u}t) + cos(uit)) 
+hna 2 (t)(2/3 2 (t) - 6 e iwt ) - hna 2 (t)b e- 



(18a) 



0. 

(18b) 



The equations in (fT7|) and (fTSj) arc ordinary differential 
equations and they are easy to be solved. However, one 
needs the theorem as given in Appendix ((B]) to simplify 
the solutions, in case they are easily modelled. 

Thus, from the equations ([17]) and (|T8j). it 
can be obtained that a(£)({2,3}) an d /?(i)({2,3}) f° r 
T>({<72, &3})[N = 1]. Since the expression for 






2000 4000 6000 8000 10000 



(a) 




(b) 



FIG. 2: The time evolution of the population of two 
states ( |2(a) I and the time evolution of Dirac potential 
in cavity quantum electrodynamics of T>(a3) for 
t = 10 5 . 



{a(*)({2,3}), /3(%2,3})} and |V ({ 2,3}) [!](*) | arc too long, 
it is not necessary to be shown in this paper. 

The potentials can be obtained by substituting those 
parameters into the equation ([7]). The plots of the Rabi 
oscillations and the potentials versus time are given in 
the figures © and respectively. 



III. CORRELATIONS OF THE ONE-DIMENSIONAL 
STATIONARY DIRAC FIELD POTENTIAL AND THE 
RABI OSCILLATIONS 

In this Section, we evaluate our results in the previous 
Section which are believed casting much new light on the 
underlying physical problem. 

The calculations in previous Section result (3i(t) 
which arc useful to obtain the transformations 
from b^b — > &j[l]6.;[l] as defined by the equa- 
tion ([5]). Furthermore, these new classical fields, 
&|[1]6 2 [1], &^[1]6 3 [1]}, are substituted into the 
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■ W(t) 



As can be seen in figure ( 1(a) I, the absence of Dirac 



0.04 



-0.02 




(a) 





(b) 



FIG. 3: The time evolution of the population of two 
states ( |3(a) ) and the time evolution of Dirac potential 
( |3(b)[ ) in cavity quantum electrodynamics of 2?(ct2) for 
* = 10 5 . 



modified Rabi frequency as expressed by the equation . 
Consequently, it influences the modified atomic inversion 
in the equation . The illustration of these schemes can 
be found on the figure ((5J) in the Appendix ([C]). 

Initially, the atom-field system has the potential, V(t), 
which the magnitude is (M7) 2 |6o| in the semiclassical 
Rabi model and (Ml) 2 (n + i) for the Jaynes-Cummings 
model. Then, it is transformed into several types of new 
potentials. The results can be divided into two parts : 
first, the transformation results the parabolic potential 
which is performed by T>{g\) ; second, the transforma- 
tions yield the harmonic oscillator potentials which can 
be achieved by I?({(T2, 0-3}). We provide in the Appendix 
(fCj) the detailed of computational methods for producing 
the graphs in this Section. 



A. The parabolic potential 

The choice of o~\ corresponds to circumvent the atomic 
population into the equilibrium condition. 



potential causes the atom inclined to occupy the ground 
state: during the potential is turned off, i.e., for t < 10, 
the first oscillations form wells meaning that initially the 
atoms occupy ground state, |<?)|n), since the wells occupy 
the negative regime. 

Furthermore, the potential is turned on after £ = 10. 
It causes that the atomic population in ground state is 
gradually vanished and followed by coherently sharing 
the atomic population between the excited- and ground- 
state due to the presence of the potential. Under the 
circumstances, the Rabi oscillations are collapsed and ex- 
tremely damped out to monotonous line. 

This phenomenon is similar to the previous results in 
Reb2 and Ref42: due to the classical drive field on res- 
onance, the population is shared coherently between the 
states \g,0) and |/,1), where |Z,1) = ^ |e ' 0) + |g ' 1>) . 

The figure (1(b)) shows that the Rabi oscillations is ~ 

—. — , since the form remembers to Yukawa's potential^, 
thus we name it by Yukawa — Rabi oscillations . This 
oscillations can be generated if the Dirac potential is 
parabolic as shown in fig (1(b) |. 



The results may be very important for the future of 
quantum information. The choice of o~\ causes the coher- 
ent share between the excited- and ground-state mean- 
ing that in the case that strong interaction of classical 
field is much greater than atomic excitation in cavity 
quantum electrodynamics, the parabolic Dirac potential 
can change the qubit basis from o^-cigenbasis into a x - 
cigenbasis. 



B. The harmonic oscillator potentials 



The figure (2(a)) shows that T>(a^) generates the se- 
quences of peaks in negative regime which means the pop- 
ulation is driven into the ground state. The sequences 
of peaks in Rabi oscillations can be created by adjust- 
ing fast oscillations of the Dirac potential in the positive 
regime as shown in the figure ( |2(b)[ ) . Under this scheme, 
the sharpness of the peak in atomic inversion is directly 
proportional to the oscillations speed of Dirac potential: 
for instance, when t ~ 2000, the extreme sharp peak is 
produced if the fast oscillations of Dirac potential occur. 
Contrarily, when t ~ 6000, the unsharp peak is caused 
because the oscillations of Dirac potential are smooth and 
slow. 

The emergence of these peaks in Rabi oscillations 
and the correlations with the Dirac potentials, may be 
relevant in the case of circuit quantum electrodynam- 
ics consisting of two superconducting qubits coupled to 
an on-chip coplanar waveguide (CPW) and explored by 
the power dependence of the heterodyne transmission as 
shown in Ref4. Bishop et al. show that the emergences 
of peaks identified with a multiphoton-transmon qubit 
transition from the ground state to an excited Jaynes- 
Cummings state can be accomplished due to the drive 
power of heterodyne transmission. In contrast with our 
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results, the transition begins to saturate as the drive in- 
creases. 

Furthermore, the choice of (02) produces Rabi oscilla- 
tions which are proportional to the oscillations of Dirac 
potential. 



It is shown in the figure (3(a) ) that initially the atom is 



in the excited state. The atomic population is gradually 
driven into the ground state along with the oscillations 
of potential gradually change from the positive- into the 
negative-regime. 

This choice demonstrates that it is possible to ensure 
the coherence of the qubits without the collapse of 
Rabi oscillations by tuning the Dirac potential in the 
oscillations form across the positive-negative regimes as 
shown in figure ( |3(b)| . This type of Rabi oscillations 
may be similar to the case of an atom which is initially 
in the excited state and field initially in a thermal stated. 



IV. PHYSICAL IMPLEMENTATION 

The proposed method suggests that the scheme is 
working under the apparatus in which the electromag- 
netic fields concurrently behave quantumly and classi- 
cally. The one of possibilities to perform this is by in- 
volving the extra device in the cavity in which it gener- 
ates the field behaving classically. Another possibility, it 
can be done by involving instrument performing external 
forces behaving as artificial classical field. 

In the recent research of photon in cavity, the atom- 
cavity system can be excited by photon transmission fol- 
lowed by the photon blockade of an optical cavity enclos- 
ing one trapped ion in the regime of strong atom-cavity 
coupling^. A single atom path is also controllable by the 
feedback of photon- by-photon^. However, it is not clear 
how photon, the quantum of the electromagnetic field, 
contributes to control the excitation and decay of the 
atomic population. 

The next constraint is the classical field contributes 
into the modified Rabi oscillations in which the magni- 
tude is ~ VWb. The possible experiment of classical field 
driving the atomic population in cavity may relate to the 
experimental research realizing the Keldysh picture^. 
Recent explorations demonstrate that controlling of elec- 
tronic motion is enabled in ultrafast laser sources in the 
mid-infrared region^. 



the Rabi oscillations transformation by involving classical 
effect of electromagnetic in the equation. This method 
shows that it is possible to control the collapse and resur- 
gence of Rabi oscillations in cavity QED under Darboux 
transformations. 

The Pauli matrices are the parameters in the 
BBGS-Darboux transformations to determine the one- 
dimensional stationary Dirac potential of electromag- 
netic field which are responsible for controlling the Rabi 
oscillations. The appropriate choice of the parameters 
may be necessary for elucidating the surprising responses 
of the oscillations due to the external perturbations as 
found in Ref^il. 

Based on the results, it is possible to propose an 
open-loop control mechanism for Rabi oscillations under 
Darboux transformations : the operator T>(pi) is as- 
sumed as the controller and the o~i is the system in- 
put. The initial system is {V, ^} and the final system 
is {V[iV], ^[AT]}. The output variables which are read 
by the sensor are the new eigenvalues of the new Hamil- 
tonian, £jv(cj). The next challenge is defining an ap- 
propriate Darboux transformations so that the theoret- 
ical explanation underlying complete closed-loop control 
mechanism, as proposed in Ref^i, can be realized. 

For further studies, it may be interesting to investigate 
the influence of A-fold Darboux transformations to the 
multi-qubits system for various types of quantum sys- 
tem in addition to Bose-Einstein condensation of exciton 
polaritonsi^, Anyons in a weakly interacting system^, 
etc. Furthermore, it is also very interesting if the Dirac 
potential is extended into n-dimensions and related to 
the various quantum systems, since the extension of the 
potential dimensions can be exploited to describe quan- 
tum transistor—. 

We expect that this paper open the possibility to in- 
volve Darboux transformations into the extensive re- 
search of quantum information. 
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Appendix A: Homotopy Perturbation Methods 

In this appendix, we provide the algorithm for obtain- 
ing exact solution of coupled nonlinear equations^. 

1. Unravel the differential equations using conven- 
tional HPM. In this step, the widely known Tay- 
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lor series is needed to expand the solution of the 
nonlinear differential equations, f(x), at the point 
x = 0. The series reads: 



(Al) 



i=0 



i=0 



{a,b,p,q € 3? and (6 ^ 0)} 



(B5) 



and sgn(b)=j^ is defined as the sign of b (to be +1 if b 
> and -1 if b <0). 

Appendix C: Computational simulation by Maple 



2. Shortened the series solution by conventional HPM. 
This step is accomplished by obtaining the embed- 
ding parameter of the solutions. 

3. Obtain the Laplace transform of the shortened se- 
ries. The well-known Laplace transform of a func- 
tion y(t) reads 



= Ms) 



y(t)e~ st dt. 



(A2) 



4. Acquire the Pade approximant of the prior step. 
The Pade approximant is a rational function to ap- 
proach the Taylor series expansion as best as pos- 
sible. It gives 



M 
~N 



M 

}{ - x) 1 + J2 ktf 

i=0 



(A3) 



where M and N are given positive integers. 
5. Obtain the inverse Laplace transform. 



Appendix B: Explicit representation for square root of a 
complex number 

Following the Refi 27 i 31 , we provide the theorem to ob- 
tain the explicit representation for square root of a com- 
plex number. 

Theorem. A complex number 



This Appendix provides the computational methods to 
generate graphs in the Section (|IIip. In this work, we use 
Maple. The following codes show the codes for modified 
atomic inversion simulation. 



> sim:= {10"ll)*eKp{-30)*sum{ ( ( (3 ) ~(n) )/{n) 1 ) *(l-( 2* (Onega) ~2)* 
sort ( Rm+n+1 ) * ( ( ( sin ( sqrt ( kappa+Rm+D+1 ) *t ) ) "2 ) / { kappa+Hm+n+1 ) 
) ) ,n«0. .1000) : 

> sol ; -subs {b [0]-1000,01,h [bar] =0 ,l,Omega=l ,omega=Pi/2, kappa= 
sqrt{2),%): 

> plot {siin{ t )+adj f t=-0, .100, numpoiiits=5, axes^f ramc,Golorfircd) } 



FIG. 4: Maple codes for the modified atomic inversion. 



The meaning of the parameters in the figure ([4]) are 
Rm is the modified Rabi oscillations, sim is the modified 
atomic inversion, adj is the adjustment for atomic inver- 
sion in case it is needed, b[0] is the amplitude of classical 
electromagnetic field, h[bar] is h, kappa is n, omega is 
10, Omega is f2, and sol is the code for setting the value 
of the constants. 

Below, the codes of {T>(<7i), T>{crz), 2?(<72)} solutions 
are shown, respectively. 



> roB:-(b[0I+({h._bar)*omcga)-l)*t+{{3*b[0]-2)*(l-{(bIbar] )* 
omega) )*si&{ omega* t ) ) / {omega* {2*b[ 0] -1) ) ; 

> pi anck ( st andard ) ; 

> reB:«(b[0]*{ (i 'bar] )*omega)-l )*t+{ (3*bl 0] -2) *{l-( {b[bar] )* 
omega) )*sin(omega*t ) ) /{omega* (2*b[0]-l) )+( ( (a "bar] )*omcga)-l) 
*b[0] *sin{omega*t)/ (2*omcga*(b[0] -1) ) ; 

> imB:=( ( (hlbar] )*Omega-l) *(2*bI0] -1 ) )/ {omega* {bl 0] -1) )+{ {3*b 
[0] -2)*{l-{h[bar] )*Omega)*cos (omega*t ) ) /(omega* (2*b[ 0] -1) )-{ < 
(h[bar] )*Omoga)-l)*bI0] *cos{omega*t)/{2*omega*{bI0] -1) )-{b[0I 
'{2)*t "{2)*omega*{ (h[bar] )*Omega-l ) )/ { { 2*bI0] -1 ) ) ; 

> absB;=sqrt{ (reB)~(2)+(imB)~2) ; 

> Kb: B siibs{bi[ 0] =10,01 , h[bar] e 0, l,omega e l, omcga e Pi ,%) ; 

> Rm: = (Hl>)~2; 



FIG. 5: Maple codes for the solution of T>(ai] 



V a + ib 



(Bl) 



can be simplified into 



by defining 



and 



where 



p + iq 



Va 2 + b 2 - a 

V2 V 



(B2) 



(B3) 



(B4) 



Ref a[3] ;t={b[0] /2 )*sin{ omega* t) } 
Imf a[3] ?={bI0] )*sin{omega*t) 

Reina[3] :=-2* {bI0] ^2)*{ { cos { omega* t) )™2) +cos{ 2*omega*t ) *{ {cos 
{omega*t ) ) ~2)*(bI0p2)-2*onsega*b[0] *cos{2*omega*t)*{sin 
{ omega* t )+cos {omega* t ) ) } 

ImibaI3] :=sin{2*omega*t )*{ { cos { omega* t ) ) ~2)*{b[0] A 2)-2* omcg a* 
b[0 ] *sin{2*omega*t )*{sin{ omega* t )+cos{omega*t ))-2*(b[0]~2)* 
sin{omega*t)*cos{omega*t) } 

Hcp[3] :=<l/sqrt(2) )*sqrt(sqrt( {Reiiia[3] ) *2+{Imlna[3] "2) )+ 
Hciba[3] ) ; 

In; [3] : = {sign{lmlna[3] )/sqrt {2) )*sgrt {sgrt { {Rcina[ 3] ) "2 + 
(Iminap] "2) )-Hclca(3] ) ; 

Rcra[3] :=Ref a [3] -Rep[3] +b[0] *sin(omega*t) ; 

Imra [ 3 ] ; =Imf a I 3 ] -Imgl 3 ] -b[ ] *cos { omega*t ) ; 

Absbb[3] :={lmra[3] ) ~2+{Rera[3] )"2 ; 

Rmt esubs (b]0] =0,1, hi bar] =0,1 ,omega=0,5 , omega=Pi , % ) j 



FIG. 6: Maple codes for the solution of 2? (03) 



Especially for 2? (172), we find the oscillation amplitudes 
are extremely huge, therefore the amplitude is scaled 
down in case it is possible to be simulated. 
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> sim; «. 9631250000- . 6806250000*008 { 3. 141592654*t)'-2+sit 

<3.111592651*t)*(l* (-cos <3.U1592651*t)"l+l*coa< 3.141592651* 

t)"2+10.*coa<3.111S92651*t)'"'6-2.*coa{3.111592651*t)'S+l*coa 

<3.111592651*tr3+9*sin<3.111592654*t)*coa<3.141592654*t)-2.* 

cos<3.141592654*t)"4*Bin<3.141592654*t)+l*(;os<3.111592651*t) 

~2*ain< 3. 111592654* t) ) -(l/2»+.3*cos( 3.141592654*t )"2- ,2*ooa 

(3.141592654*t)''4+.l*cos(3.141592654*t)"2*aill{3.141592654*t> 

+.l*coa(3.111592651*t)*3-.6*ain(3.111592651*t)-.6*coa 

<3.141592654*t))"<l/2)-.6*aill<3.141592654*t)*cos<3.141592651* 

t)+.2*(4-4*coa<3.141592651*t) -l+5*co8{3.141592654*t)*2+10.* 

cos{3.141592654*t)'6-2.*coa(3.141592654*t)"5+l*cos 

{3.141592654*t)~3+9*sill(3.141592654*t)*coa(3.141592654*t)-2.* 

cos{3.141592654*t)'4*aiD<3,141592654*t)+l*cos<3,141592654*t) 

~2*ain(3.141592654*t) )-<l/2»-.3*{l*<4-4*cos<3.111592651*t) 

-4+5*cos<3.141592654*t)"2+10.*cos(3.111592654*t)"6-2.*cos 

<3.141592654*tr5+l*cos<3.141592654*t)~3+9*ain<3.141592654*t) 

*cos(3,141592651*t)-2,*cos<3,141592651*t) "4*sin< 3. 141592654* 

t)+l*cos(3.111592654*t)-2*alll(3,141592654*t))*(l/2)+,3*ooa 

< 3, 141592654*t)' 2-, 2*coa< 3, 141592654*t)~4+,l*cos (3,141592651* 

t)'-2*aln(3.141592654*t)+.l*cos(3.141592654*t)~3-.6*sin 

<3.141592654*t»-.6*coa<3.141592654*t))~<l/2)*<:o8(3.141592651* 

t)-.5*aill(3.141592651*t)*(10*(4-4*coa{3.1415926S4*t)~4+5*cos 

<3.141592654*t)"2+10.*coa(3.141592654-t)~6-2.*<:oB 

(3.141592654*t)--5+l*ooa(3.141592654*t>~3+9*aill(3.1415926S4*t) 

*coa(3.141592651*t)-2.*ooa(3.141S92651*t)'l*aln(3.141592654* 

t)+l*coa<3.141592654*t)"2*8in<3.141592654*t) ) ~<l/2)-coa 

< 3. 141592654*t)"2+.2*<:os< 3. 141592654*t) "4-. l*cos (3.141592651* 

t) "2*sill(3.141592651*t)-.l*cos(3.141592654*t) ~34-.6*sill 

<3.141592654*t)-l-.6*cos<3.141592654*t))~(l/2); 



FIG. 7: Maple codes for the modified atomic inversion 
in the case of T>(a2)- 




FIG. 8: The scheme for obtaining the correlations 
between atomic inversion and Dirac potentials. 



